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We study the kinetic and chemical equilibration in "infinite" parton matter within the parton- 
hadron-string dynamics off-shell transport approach, which is based on a dynamical quasiparticle 
model (DQPM) for partons matched to reproduce lattice QCD results — including the partonic 
equation of state — in thermodynamic equilibrium. The "infinite" parton matter is simulated by 
a system of quarks and gluons within a cubic box with periodic boundary conditions, at various 
energy densities, initialized out of kinetic and chemical equilibrium. We investigate the approach 
of the system to equilibrium and the time scales for the equilibration of different observables. We, 
furthermore, study particle distributions in the strongly interacting quark-gluon plasma (sQGP) 
including partonic spectral functions, momentum distributions, abundances of the different parton 
species, and their fluctuations (scaled variance, skewness, and kurtosis) in equilibrium. We also 
compare the results of the microscopic calculations with the ansatz of the DQPM. It is found that 
the results of the transport calculations are in equilibrium well matched by the DQPM for quarks 
and antiquarks, while the gluon spectral function shows a slightly different shape due to the mass 
dependence of the gluon width generated by the explicit interactions of partons. The time scales for 
the relaxation of fluctuation observables are found to be shorter than those for the average values. 
Furthermore, in the local subsystem, a strong change of the fluctuation observables with the size 
of the local volume is observed. These fluctuations no longer correspond to those of the full system 
and are reduced to Poissonian distributions when the volume of the local subsystem becomes small. 
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I. INTRODUCTION 



Nucleus-nucleus collisions at ultrarelativistic energies 
are studied experimentally and theoretically to obtain 
information about the properties of hadrons at high den- 
sity and/or temperature as well as about the phase tran- 
sition to a new state of matter, the quark-gluon plasma 
(QGP). Whereas the early "big bang" of the universe 
most likely evolved through steps of kinetic and chem- 
ical equilibrium, the laboratory "tiny bangs" proceed 
through phase-space configurations that initially are far 
from an equilibrium phase and then evolve by fast expan- 
sion. On the other hand, many observables from strongly 
interacting systems are dominated by many-body phase 
space such that spectra and abundances look "thermal." 
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It is thus tempting to characterize the experimental ob- 
servables by global thermodynamical quantities such as 
"temperature," chemical potentials or entropy [il-Q- We 
note that the use of macroscopic models such as hydrody- 
namics |9l-[l2| employs as a basic assumption the concept 
of local thermal and chemical equilibrium in the infinite- 
volume limit, although by introducing different chemical 
potentials one may treat chemical off-equilibrium also 
in hydrodynamics. The crucial question, however, of 
how and on what time scales thermodynamic equilibrium 
can be achieved is presently a matter of debate. Thus 
nonequilibrium approaches have been used in the past to 
address the problem of time scales associated with global 
or local equilibration |13i-i20'|. Another question is the 
influence of finite-size effects on fluctuation observables 
and the possibility of relating experimental observations 
in relativistic heavy-ion collisions to the theoretical pre- 
dictions obtained in the thermodynamic limit. Therefore, 
a thorough microscopic study of the questions of ther- 
malization and equilibration of confined and deconfined 



matter within a nonequilibrium transport approach, in- 
corporating both hadronic and partonic degrees of free- 
dom and the dynamic phase transition, appears timely. 
The paper is organized as fohows. In Sec. II we pro- 
vide a brief reminder of the ofF-sheh dynamics and the 
ingredients of the transport approach. We then present 
in Sec. Ill the actual results on the chemical equilibra- 
tion of the partonic matter in parton-hadron-string dy- 
namics (PHSD). In Sec. IV we investigate the properties 
of the partonic matter in chemical and kinetic equilib- 
rium and compare the particle properties in equilibrium 
with the dynamical quasiparticle model (DQPM), which 
has been developed to describe the thermodynamics of 
lattice QCD. In Sec. V we study (within the dynamical 
approach) the parton properties at finite quark chemical 
potential fiq, while in Sec. VI higher moments of parton 
distributions and the equilibration of fluctuation observ- 
ables as well as the size of fluctuations in equilibrium are 
investigated. We then show in Sec. VII the time scales for 
the relaxation of fluctuation observables in comparison to 
the time scales for the equilibration of the average values 
of the observables. Finally, a summary and conclusions 
are given in Sec. VIII. 



II. THE PARTON-HADRON-STRING 
DYNAMICS TRANSPORT APPROACH 

In this work we study the kinetic and chemical equili- 
bration in "infinite" parton matter within the Parton- 
Hadron-String Dynamics (PHSD) transport approach 
[21I |22| . which is based on generalized transport equa- 
tions on the basis of the off-shell Kadanoff-Baym equa- 
tions (23, UM fo'^ Green's functions in phase-space rep- 
resentation (in the first order gradient expansion, be- 
yond the quasiparticle approximation). The approach 
consistently describes the full evolution of a relativistic 
heavy-ion collision from the initial hard scatterings and 
string formation through the dynamical deconfinement 
phase transition to the strongly interacting quark-gluon 
plasma (sQGP) as well as hadronization and the subse- 
quent interactions in the expanding hadronic phase. In 
the hadronic sector PHSD is equivalent to the Hadron- 
String-Dynamics (HSD) transport approach [2^, [2g that 
has been used for the description of pA and AA colli- 
sions from GSI heavy ion synchrotron (SIS) to Relativis- 
tic Heavy Ion Collider (RHIC) energies in the past. 

In particular, PHSD incorporates off-shell dynamics 
for partons and hadrons. In the off-shell transport 
description, the hadron and parton spectral functions 
change dynamically during the propagation through the 
medium and — in case of hadrons — evolve toward the on- 
shell spectral function in vacuum if the system expands in 
the course of the heavy- ion collisions. As demonstrated 
in [27|,[2g| the off-shell dynamics is important for hadronic 
resonances with a rather long lifetime in vacuum but 
strongly decreasing lifetime in the nuclear medium (es- 
pecially u and (j) rnesons) and also proves vital for the 



correct description of dilepton decays of p mesons with 
masses close to the two-pion decay threshold. 



A. Off-shell transport 

Let us recall the off-shell transport equations (see [23| 
for details) . One starts with a first-order gradient expan- 
sion of the Wigner-transformed Kadanoff-Baym equation 
and arrives at the generalized transport equation |2J, [23] 

2p^9yG>< - {ReE^,iG><}-{iS><,ReG^} 



{M, iG><} 



it<iG> -it>iG< , (1) 



where the curly brackets denote the relativistic general- 
ization of the Poisson bracket 

{F,G} = dlF{p,x)d^,G{p,x) - d^F{p,x)dPGip,x) . 

One additionally obtains a generalized mass-shell equa- 
tion 
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= i t>< ReG« + \{iJ:>,iG<} - \{iI]<,iG>} (2) 

with the mass function M = p^ — m? — ReS''^. In Eqs. 
(HI) and ([2]) the Green's functions G^^ stand for the ex- 
pectation values of the quantum fields (denoted here by 



i G^(x,y) = 7/(<i>"f(y)$(a;)) 
zG>(a;,y) = (a>(x)$t(y)), 



(3) 



with ?7 = 1 for bosons and 77 = —1 for fermions, while the 
self-energies S(a;, y) are given by the functional derivative 
of J- with respect to the full propagator G: 



2i- 
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(4) 



In (JH) the functional J- is the sum of all closed two- 
particle-irreducible (2PI) diagrams built up by full prop- 
agators G. 

The retarded and advanced Green's functions G^ and 
G^ are given as 

G''{XI,X2) = e(ii - t2) [G>{XI,X2) ~ G<{XI,X2)\ , (5) 
G^{xi,X2)^-Q{t2-t^)[G>{xi,X2)-G<{xi,X2)\.{&) 

These Green's functions contain exclusively spectral but 
no statistical information of the system. Their time evo- 
lution is determined by Dyson- Schwinger equations (cf. 
Ref. 27|). Retarded and advanced self-energies S^ and 
T,^ are defined in analogy to Eqs. ^ and ([5]). 

In the transport equation ([T]) one recognizes on the 
left-hand side the drift term p^d^ iG^"^, as well as the 



Vlasov term with the real part of the retarded self-energy 
ReS-^. On the other hand the right-hand side represents 
the collision term with its typical "gain and loss" struc- 
ture. Thus interaction between the degrees of freedom 
is thus incorporated into the mean fields and collisions 
as in the Vlasov-Boltzmann "standard" transport ap- 
proach w]. In contrast, in the off-shell transport there 
is an additional term —{i'S^^jReG^}, which is denoted 
as the back-flow term and is responsible for the proper 
off-shell propagation. It vanishes in the on-shell quasi- 
particle limit. Note, however, that the self-energies S 
fully determine the dynamics of the Green's functions 
for given initial conditions. 

We, further on, represent Eqs. (H]) and ([2]) in terms of 
real quantities by the decomposition of the retarded and 
advanced Green's functions and self-energies as 

qB/A ^ ^^qR ^ ■ ^^qB. ^ ^^qB. zp iA/2 , (7) 

A^T^lmG"/^, (8) 

S«M = ReE« ± ilmE^ = ReS« T «f/2 , (9) 
f = =F2Imf;«/^ . (10) 

We note that in Wigner space the real parts of the re- 
tarded and advanced Green's functions and self-energies 
are equal, while the imaginary parts have opposite sign 
and are proportional to the spectral function A and to 
the width F, respectively. 

With the redefinitions ([7l)- p^ one obtains two alge- 
braic relations for the spectral function A and the real 
part of the retarded Green's function, ReG^, in terms of 
the width F and the real part of the retarded self-energy, 
ReE^, as [13, ^ 

[pI -p^-m^- ReE^] ReG^ = 1 + | f A , (11) 
[pI - p2 - m2 - ReE^] A = f RcG^ . (12) 

Note that all terms with first-order gradients have disap- 
peared in (TTlT) and p2|) . A first consequence of ([12]) is a 
direct relation between the real and the imaginary parts 
of the retarded, advanced Green's function, which reads 
(for f ^ 0) 



ReG'^ = 



pI 



p2 - m^ - ReE^ 



(13) 



Inserting (|13p in (jlip we end up with the following result 
for the spectral function and the real part of the retarded 
Green's function: 



A 



[pI 



RcE' 



r2/4 



ReG^ = 



[pI 



ReE' 



(14) 
(15) 



The spectral function ([14]) shows a typical Breit- Wigner 
shape with energy- and momentum-dependent self- 
energy terms. Although the above equations are purely 



[pI - P^ - m^ - ReE^] -t- rV4 



algebraic solutions and contain no derivative terms, they 
are valid up to the first order in the gradients. 

In addition, subtraction of the real parts and adding up 
the imaginary parts lead to the time evolution equations 

pt'd'^^A^^ {Ret^,A } + ^ {f,ReG^} , (16) 
p'^a^^ReG^ = I {ReS«,ReG^} - | {f, A} . (17) 

When inserting ([H]) and ([15]) we find that these first- 
order time-evolution equations are solved by the algebraic 
expressions. Accordingly, the time evolution of the sys- 
tem is fully defined by ReE^ and the width F in ([1]). 

We recall that the off-shell transport equation (1) can 
be solved explicitly by employing a generalized test- 
particle ansatz for the real quantity JG^(x,p). For the 
explicit equations of motion for these test particles we 
refer the reader to Ref . i29:j . 



B. Explicit equations for fermions 

In case of fermions — such as baryons or quarks — 
the self-energy ReE^ is separated into different Lorentz 
structures of scalar and vector type: 

ReE^/m,, = U^{x,p) + j^Ui:{x,p) (18) 

for each fermion species h. The mass function for 
fermions is then 



Mh{p,x)^nl 



n^ 



*2 



(19) 



with the effective mass and four-momentum given by 



ml{x,p) = nih + U^{x,p) , 
U^'{x,p)^p^'-Ui:{x,p) , 



(20) 
(21) 



where rrih stands for the bare (vacuum) mass. After in- 
serting (|19p into the generalized transport equation ([TJ, 
the covariant off-shell transport theory emerges. It is for- 
mally written as a coupled set of transport equations for 
the phase-space distributions Nh{x,p) [x = (i, r), p ~ 
(w, p)] of fermion h with a spectral function Ah{x,p) [us- 
ing iG^{x,p) = Nh{x,p)Ah\x,p)], i.e., 

(n^ - n,a^[/,r - midiu^) d^N,,{x,p)A^{x,p) 
+ (n,a^f/,r + midiui) d^Nj,ix,p)Ahix,p) 

-{it<,ReG^} 
= (27r)4^ ir2tr3ir4[rtT]i2^34<5^(n+n2-n3-n4) 

X Ah{x,p)Ah:,ix,P2)Ah3{x,P3)Ahi{x,P4) 
X [Nh^{x,P3)Nh^{x,P4)fh{x,p)fh2{x,P2) 

-Nh{x,p)Ni,2{x,p2)fh3{^,P3)fin{x,P4)] (22) 



with 



fhix,p) = 1- Nh{x,p) 



and 



trn = 



(27r)4 



Here d^ = (9t,Vr) and df^ 



(a^,Vp) ifi = 0,1,2,3). 
The factor |T''^T| stands for the in-medium transition ma- 
trix element (squared) for the binary reaction 1 + 2 — > 
3-1-4, which has to be known also off the mass shell. The 
back-flow term in (1221) is given by 



-{iS<,ReG^} 
-Of. 



Mh{x,p) 



P\Mh{x,pY+Th{x,pY/A) ^ 
Mh{x,p) 



dUNh{x,p)rh{x,p)] 



(23) 



As pointed out before this expression stands for the off- 
shell evolution, which vanishes in the on-shell limit or 
when the spectral function Ah{x,p) does not change its 
shape during the propagation through the medium, i.e., 
for VrT{x,p) = and VpT{x,p) = 0. We recall that the 
transport equation ([2^ has been the basis for the off-shell 
HSD transport approach for the baryon and antibaryon 
dynamics. 

In order to specify the dynamics of partons one has 
to evaluate/specify the related self-energies for quarks 
and antiquarks as well as gluons that enter the spectral 
functions ([T4l) and retarded Green's functions (fT5|) . This 
task has been carried out within a dynamical quasiparti- 
cle model. 



C. The dynamical quasiparticle model 

The basis of the partonic phase description is the dy- 
namical quasiparticle model [3lL |32], which has been 
matched to reproduce lattice QCD results (IQCD) — 
including the partonic equation of state — in thermody- 
namic equilibrium. The DQPM allows for a simple 
and transparent interpretation of thermodynamic quan- 
tities as well as correlators — measured on the lattice — by 
means of effective strongly interacting partonic quasipar- 
ticles with broad spectral functions. The essential quan- 
tities in the DQPM are "resummed" single-particle prop- 
agators Gq, Gq, and Gg. We stress that a nonvanishing 
width r in the partonic spectral function is the main dif- 
ference between the DQPM and conventional quasiparti- 
cle models |33i] . Its influence on the collision dynamics is 
essentially seen in the correlation functions; e.g., in the 
stationary limit, the correlation involving the off-diagonal 
elements of the energy-momentum tensor T'^' define the 
shear viscosity rj of the medium [SJ] • Here a sizable width 
is mandatory to obtain a small ratio of the shear viscosity 
to entropy density, rj/s [35| , which results in a roughly hy- 
drodynamical evolution of the partonic system in PHSD 
[22i] . The finite width leads to two-particle correlations. 



which are taken into account in PHSD by means of the 
generalized off-shell transport equations (cf. Sec. II. A) 
that go beyond the mean-field or Boltzmann approxima- 
tions. 

In the scope of the DQPM the running coupling con- 
stant (squared) for T > Tc is approximated by 



9\T/T,) 



487r2 



(117V, - 27V/) ln[A2(T/T, - T,/T,)^] ' 



(24) 



where the parameters A = 2.42 and Tg/Tc — 0.56 have 
been extracted from a fit to the lattice QCD equation 
of state as described in Refs. [M El- In ^, N^ = 3 
stands for the number of colors, T, is the critical temper- 
ature (=158 MeV), while 7V/(=3) denotes the number of 
flavors. 

In the asymptotic high-momentum (high-temperature) 
regime, the functional form of the parton quasiparticle 
mass is chosen to coincide with that of the perturbative 
thermal mass, i.e., for gluons 

7V/,^(T) = ^(^(7V. + l7V,)r^-f^^|j , (25) 



and for quarks (antiquarks) 



a' T^ ' "' 



2 I ' 



(26) 



but with the coupling given in ([24l) . The effective quarks, 
antiquarks, and gluons in the DQPM have finite widths, 
which for fiq = are approximated by 



r,(,-)(T) 



1 g^T 

3 Stt 

ITV^ 

3 27V, Stt 



ig'T 



2c 
In 



1 



(27) 
(28) 



where c = 14.4 (from Refs. [3J]) is related to a magnetic 
cutoff. Note that for /x^ = the DQPM gives 



M, 



q{q) 



3 ^' 



^n(n) = -Po 



«(«) 



9 



(29) 



From the expressions ([M|) - (l^^ . one can see that at high 
temperature, T — > oo, the masses and the interaction 
strength of the quasiparticles in the DQPM are approach- 
ing the one-loop perturbative QCD results. However, the 
one-loop functional form is not the relevant description 
at temperatures close to T, or even below. The transition 
region (approximately 0.9 Tc < T < 1.1 T,) is dominated 
by nonperturbative phenomena. Therefore, we imple- 
ment in PHSD for the transitional values of T/Tc (the 
region 0.9rc < T < LIT,) functional forms for Mg^g and 
Tqg, which are growing softer with decreasing T/Tc as 
compared to the perturbative logarithmic divergence in 
([M)l . The actual values of Mq^g and Fg^^ have been shown 
as functions of temperature as well as the scalar parton 



density in Ref. |4(]||. A comparison to the IQCD interac- 
tion measure [33| has been presented also in Ref. 43] (cf- 
Figs. 1, 2, and 4). 

We note in passing that the smooth parametrizations 
for Mq^g and Vq^g at T close to Tc nicely reproduce the 
recent IQCD calculations from the Wuppertal-Budapest 
group 39]. As a consequence, we obtain not only a 
quantitatively good description of the phase transition re- 
gion but also a smooth "interpolation" from the hadron- 
dominated systems to those with dominant partonic de- 
grees of freedom. 

With the parton masses and widths fixed by ((24)) - (p9)) 
the spectral functions can be written [in alternative form 
to ([Tl)l as 



Aj = pj(w,p)= — 



EA{io-E,f+n {uj + E,f+T] 



AujT. 



(tj2 - p2 _ M]Y + 4r2tj2 ' 



(30) 



separately for quarks, antiquarks, and gluons {j = 
q,q,g), with the notation -B|(p^) = p^ + Mf - T^. We 
may identify (cf. Sec. II. A) 



-ReEf = Mf, r,- = 2ujTj 



(31) 



The spectral function (|30p is antisymmetric in w and nor- 
malized as 



duj f dio 

— ujpj{uj,p) = / ^ 2wpj(w,p) ^ 1 . (32) 



The parameters Tj and Mj from the DQPM have been 
defined above in the Eqs. (25)-(29). Note that the 
DQPM assumes Tj — const(a;); we will discuss the con- 
sequences of this approximation in Sec. IV. Also, the 
decomposition of the total width Tj into the collisional 
width (due to elastic and inelastic collisions) and the de- 
cay width is not addressed in the DQPM. Therefore, we 
dedicate the next section to this question and to a brief 
description of the microscopic implementation of the de- 
confined phase of QCD within the PHSD. 



D. Reaction rates and effective cross sections 

In this section we present the effective cross sections 
for each of the various partonic channels as a function 
of energy density e; these cross sections determine the 
partial widths of the dynamical quasiparticles as well as 
the various interaction rates. This analysis is important, 
because, although the DQPM provides the basis of the 
description of the strongly interacting quark-gluon sys- 
tem in PHSD in equilibrium, the dynamical transport 
approach (i.e., PHSD) goes beyond the DQPM in sim- 
ulating hadronic and partonic systems also out of equi- 
librium. For the microscopic transport calculations, the 



partial widths of the microscopic scattering and decay 
channels have to be known, while the DQPM provides 
only the total widths of the dynamical quasiparticles 
that have been fixed by lattice QCD calculations as de- 
scribed in Sec. II. C and in more detail in Refs. p3ll.l33|. 
Furthermore, the explicit shape of the partonic spectral 
functions — taken as Lorentzians in the DQPM (PU)) — will 
depend on the decomposition of the interaction into par- 
ticular channels within the coupled-channel dynamics of 
PHSD. 

In order to fix the partial cross sections for the interac- 
tions between the dynamical quarks and gluons (as func- 
tions of energy density e) we perform PHSD calculations 
in a cubic finite box with periodic boundary conditions — 
simulating "infinite" hadronic or partonic matter. In 
this particular case the derivatives of the retarded self- 
energies with respect to space vanish in (|22p such that 
we essentially deal with the parton dynamics due to the 
collision terms in ([2^ . 

The following (quasi)elastic interactions among 
quarks, antiquarks, and gluons (g, q, g) are implemented 
in PHSD: 



q{mi)+q{m2) -^ q{mi,) + q{mi) , 


(33) 


q + q -^ q + q, 


(34) 


q + q -^ q + q, 


(35) 


g + q -^ g + q, 


(36) 


g + q -^ g + q, 


(37) 


9 + 9 -^ 9 + 9- 


(38) 



The (quasi)elastic processes (l33t - (|38|) play a crucial role 
for the thermalization in PHSD due to the possibility to 
change the masses of interacting partons in the final state 
as shown in Eq. ([55t . 

The flavor exchange of partons is possible only within 
the inelastic interactions in PHSD, which are: 



9 ^ q + q, 

9 ^ 9 + 9, 
g + g O q + q. 



(39) 
(40) 
(41) 



The inelastic interactions (|39p - ((lT|) are the basic pro- 
cesses for the chemical equilibration in PHSD; however, 
the inelastic processes [(HOI and (HD)] are strongly sup- 
pressed (<1%) kinematically in PHSD due to the large 
masses of gluons. 

We recall that for binary channels we have explicit for- 
mulas for the partial widths, e.g. [from the collision term 
in (Ell)], 



r-'-*-(pi) = ^ tr2ir3tr4|rtT|2^2^; 



3+4 



2,3,4 



X ^^2 {p2)Ah3 {p:i)Ah^ {pi)Nh., (P2)//i3 (P3)//i4 iPi) 

x{2nfS\Pi+P2-P3-P4) , (42) 

where hi is an index, which can be equal to "gi" , "g^" or 
"(7i," where i = 1,2,3. Since we study partons at high 



temperature the fermion blocking terms can be neglected, 
i.e., approximated by / = 1, and one ends up with 



elastic 



iPl) 



E 

2,3,4 



tr2 \T^T\l 



+2^3+4 



xAh^{p2)Nh^{p2)R2ipi+P2;M3,Mi) , (43) 

where the four- momenta of particle 4 are fixed by energy- 
momentum conservation and R2 stands for the two-body 
phase-space integral (cf. [41|). We recall that the squared 
matrix element times the two-body phase-space integral 
defines a binary cross section a times a kinematic factor, 
i.e.. 



E\T'T\l 



+2-!-3+4 



R2iP3 + Pi) = 4.EiE2Vrel<J, (44) 



3,4 



with the relativistic relative velocity for initial invariant 
energy squared, s, given by 



Vrei = yis- Ml - Ml) - 4Mf M|/(2^i^2) . (45) 

In (|44p X^s 4 stands for a summation over discrete final 
channels. 

If the cross section a is essentially independent of the 
momenta, which should hold for low-energy binary scat- 
tering, we may write (1431) as 



-^elastic 



(pi) = {vi2a)N2 



(46) 



which corresponds to the Boltzmann limit relating the 
collision rate to the average velocity between the collid- 
ing partners (in the center-of-mass frame) and the cross 
section for scattering as well as the density N2 (summed 
over the discrete quantum numbers of particle 2). We 
employ these relations in determining the effective elastic 
cross sections between partons in the PHSD. Note that 
the total number of collisions between particles of type 1 
and 2 are obtained from (^S)) (in our case) by multipli- 
cation with the volume V and the particle density TVi, 



dt 



V{vi2cr)NiN2 



(47) 



Both the number of collisions between the individual par- 
ticle species as well as their densities are easily accessible 
in the transport approach. 

The cross section for gluon formation from fiavor- 
neutral q + q interactions in the color octet channel is 
calculated by the resonant cross section at invariant en- 
ergy squared, s = {pq +Pq)'^, 



(Tqq-,g{s,e,Mq,Mq) = - 



47rsr2(£) 



1 



[s~M^{e)Y + srl{s)Pr'ei' 



with 



p2 [s-{Mq + Ml] [s~iMq-Mqf] 



rel 



As 



(48) 
(49) 
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FIG. 1: (Color online) The energy density dependence of the 
coefficient aqq extracted from the PHSD simulations in the 
box (blue dots) and corresponding fit (red line). 



while the factor 2/4 corresponds to the ratio of final to 
initial spins (assuming two transverse degrees of freedom 
for the gluon in line with the DQPM). Note that formula 
(l48t provides an off-shell cross section which depends on 
the four-momenta of the incoming quark and antiquark 
as well as on the spectral properties of the gluon. We 
recall that in the actual simulation the quark and anti- 
quark masses are distributed according to the spectral 
function pO|) and their three-momenta vary in a broad 
range roughly in line with thermal Boltzmann distribu- 
tions. 

We point out that the iteration of the coupled equa- 
tions has been performed with the additional boundary 
conditions 



'^gq(qg) ~ g'^'ssl^)' 



(^qq{£)<^gg{£) (50) 



as suggested by lattice QCD, which roughly follows a 
scaling with the color Casimir operators. This is also 
refiected in the DQPM ansatz (P^ . We mention that 
this scaling might be violated and require a further inde- 
pendent parameter, which, however, presently cannot be 
fixed appropriately by IQCD calculations. The function 
aqq{e) has to be determined by the iteration procedure 
until self-consistency has been reached for each value of 
energy density e. Note that for /ig = we have identical 
phase-space distributions for quarks and antiquarks and 
also identical interaction rates, which simplifies substan- 
tially the iteration process. Additionally, we assume for 
the present study that the elastic scattering process is 
isotropic. 

The numerical results of the self-consistent determina- 
tion of the cross sections and widths can be parametrized 
in the following form (with the cross sections given in 
square femtometers) : 

(Tgg{s) w 7.6e-"/"-8+106.2e^"/0-2+1.7e-^/3.7^o 3^ (5^) 
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FIG. 2: (Color online) The gluon-gluon (solid blue line), 
quark(antiquark)-quark(antiquark) (dashed red line), and 
quark(antiquark)-gluon or gluon-quark(antiquark) (dash- 
dotted green line) elastic cross sections as functions of the 
energy density. 



where e is given in units of GeV/fm^. The solution of 
the coupled equations then give the coefficient 



aggie) wO.3 e"=/2.6_^Q4 



(52) 



This fit is shown in comparison to the numerical results 
of the iteration in Fig. [T] Accordingly, the expressions 
for the partonic elastic scatterings may be parametrized 



'^ggigg) 



,{e), agg^{0.3e~'/^-'' + 0A)agg{e). (53) 



In Fig. [2] we display the resulting gluon-gluon (solid blue 
line), quark-quark (dashed red line), and quark- gluon 
(dash-dotted green line) elastic cross sections as func- 
tions of the energy density. Note that these cross sec- 
tions are moderate at high energy density and typically 
in the order of 2-3 mb but become large close to the crit- 
ical energy density. This behavior basically reflects the 
infrared enhancement of the strong coupling p4)) around 
Tc and implies that partons "see each other" at distances 
of about 1 fm (and even more) in the vicinity of the phase 
transition. The physics interpretation is that color sin- 
glet qq pairs form "rotating strings" whereas qq or (qq) 
pairs form resonant (and colored) diquark (antidiquark) 
states that may fuse with another quark (or antiquark) 
to form baryonic resonances. 

Although the cross sections ([5^ have been extracted 
for fig = in thermal equilibrium we may adopt the same 
cross sections also out of equilibrium and for /i^ ^ in 
the PHSD transport approach. This appears legitimate 
for phase-space configurations slightly out of equilibrium 
as well as for moderate fig. 



E. Dynamical hadronization 

In the present manuscript we essentially consider 
systems in the partonic phase where the dynamical 
hadronization plays no substantial role. However, we 
describe here in short the implementation of the tran- 
sition from the partonic to hadronic degrees of freedom 
(hadronization) and vice versa (deconfinement) in PHSD. 
Hadronization is described in PHSD by covariant tran- 
sition rates for the fusion of quark-antiquark pairs to 
mesonic resonances or three quarks (antiquarks) to bary- 
onic states [22|, e.g., for g-)-q fusion to a meson m of four- 
momentum p = (lu, p) at space-time point x = (t, x): 



dNrn{x,p) 

d'^xd'^p 



--TrgTrgS\p-Pg-Pg)S'l'''^+'''' 



XUJgPg{Pg)ujgPg{Pg)\Vgg\'^Wm ( Xg " Xg, "^ '^ 



xNg{xg, Pg)Ng{xg, Pg)6(JieivOT, color 



(54) 



In (I54p we have introduced the shorthand operator nota- 
tion 



Tr, 



E 



d Xi 



^P, 



{2^Y 



(55) 



where ^ denotes a summation over discrete quantum 
numbers (spin, flavor, and color); Nj{x,p) is the phase- 
space density of parton j at space-time position x and 
four-momentum p. In (1541) (5(flavor, color) stands sym- 
bolically for the conservation of flavor quantum numbers 
as well as color neutrality of the formed hadron m, which 
can be viewed as a color dipole or "pre-hadron." Further- 
more, Vqq{pp) is the effective quark-antiquark interaction 
from the DQPM (displayed in Fig. 10 of Ref. i31lj) as a 
function of the local parton {q + q + g) density pp (or en- 
ergy density). Furthermore, Wm{x^p) is the dimension- 
less phase-space distribution of the formed pre-hadron, 



Wm{S.,P() = exp — 2 exp 



262 



2V 



v\ 



{Mg 



Mg? 



(56) 



with ^ = xi — X2 = Xq — Xq and pj = (pi — P2)/2 = 
{j)q — pq)/2. The width parameter b has been flxed by 
\J {r'^) = b = 0.66 fm (in the rest frame), which corre- 
sponds to an average rms radius of mesons. We note 
that the expression ([5S|) corresponds to the limit of in- 
dependent harmonic oscillator states and that the final 
hadron-formation rates are approximately independent 
of the parameter b within reasonable variations. By con- 
struction the quantity (j56p is Lorentz invariant; in the 
limit of instantaneous "hadron formation," i.e., ^^ = 0, 
it provides a Gaussian dropping in the relative distance 
squared, (ri — r2)2. The four- momentum dependence 
reads explicitly (except for a factor 1/2) 

(El - E2f - (pi - P2)' - (Afi - Nhf < (57) 
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FIG. 3: (Color online) PHSD calculations for the system 
initialized by quarks and gluons a.t ^q — and e — 0.35 
GeV/fm''. The numbers of partons (solid red line) and 
hadrons (dashed blue line) are shown as functions of time. 



and leads to a negative argument of the second expo- 
nential in ()56p favoring the fusion of partons with low 
relative momenta pq — pq — pi ~ p2- 

Note that due to the off-shell nature of both partons 
and hadrons, the hadronization process obeys all conser- 
vation laws (i.e., four-momentum conservation and flavor 
current conservation) in each event, the detailed balance 
relations, and the increase in the total entropy S in case 
of a rapidly expanding system. The physics behind (I54p is 
that the inverse reaction, i.e., the dissolution of hadronic 
states to quark- ant iquark pairs (in case of mesons), at 
low energy density is inhibited by the huge masses of 
the partonic quasiparticles according to the DQPM. Vice 
versa, the resonant q-q pairs have a large phase space to 
decay to several 0~ octet mesons. We recall that the 
transition matrix element becomes huge below the criti- 
cal energy density [32j. For further details on the PHSD 
ofF-shell transport approach and hadronization we refer 
the reader to Refs. |2l|, [H, ^28, 29, 40]. 

If the system is initialized by an ensemble of partons, 
but the energy density in the system is below the crit- 
ical energy density {sc ~ 0.5 GeV/fm'^), the evolution 
proceeds through the dynamical phase transition (as de- 
scribed in Sec. IIIEp and ends up in an ensemble of 
hadrons. In Fig.|3]we show the results of the PHSD calcu- 
lations for the system initialized by quarks and gluons at 
fiq = and e = 0.35 GeV/fm'^. The numbers of partons 
(solid red line) and hadrons (dashed blue line) are shown 
as functions of time. We observe that the transition from 
partonic to hadronic degrees of freedom is complete after 
about 9 fm/c. A small nonvanishing fraction of partons 
remains due to local fluctuations of energy density from 
cell to cell. The equilibration of hadron-dominated mat- 
ter is an interesting topic. However, we concentrate in the 
present work on the properties of parton-dominated mat- 



ter, since we are primarily interested in the time scales 
for kinetic and chemical equilibration in the sQGP within 
the PHSD. Thus we will study the systems at energy den- 
sities higher than the critical one for the remainder of this 
work. 



III. CHEMICAL AND THERMAL 
EQUILIBRATION 

Before we proceed to the actual results on the chemical 
equilibration and kinetic thermalization of the partonic 
matter in PHSD, let us note that the PHSD transport 
approach has been tested in comparison to various data 
from relativistic heavy-ion collisions and has led to a fair 
description of particle production 21 1, elliptic flow JSTJ, 
and dilepton production [3^ both at Super Proton Syn- 
chrotron (SPS) and top RHIC energies. In particular, 
the comparison of PHSD calculations to the data of the 
NA60, PHENIX, and STAR Collaborations in Ref. H 
has shown that the partonic dilepton production channels 
should be visible in the intermediate-mass region (from 
1 to 3 GeV). The partonic contribution to the dilepton 
radiation appears to be exponential in mass from 1 to 2.5 
GeV so that an interpretation in terms of thermal radi- 
ation from the sQGP might appear appropriate. How- 
ever, such an interpretation is subject to the question 
of whether or not the PHSD dynamics shows that ki- 
netic equilibrium is achieved on the partonic level within 
the characteristic lifetime of the partonic system in these 
collisions. We will address this question in the present 
section. 

As mentioned above, we simulate the "infinite" mat- 
ter within a cubic box with periodic boundary conditions 
at various values for the quark density (or chemical po- 
tential) and energy density. The size of the box is fixed 
to 9^ fm^ for most of the following calculations. How- 
ever, we will study also larger box sizes in order to deter- 
mine whether the thermodynamic limit is approximately 
reached, in particular when addressing the fluctuation 
measures. The initialization is done by populating the 
box with light (m and d) and strange (s) quarks, an- 
tiquarks, and gluons. The system is initialized out of 
equilibrium and approaches kinetic and chemical equilib- 
rium during its evolution by PHSD. We are not interested 
here in very far nonequilibrium configurations, such as, 
for example, the result of the initial hard scatterings in 
a heavy-ion collision. Instead, we study here configura- 
tions which are reasonably close to equilibrium, because 
in this case the approach to equilibrium will have uni- 
versal characteristics that will not depend on the precise 
choice of the initial state. We will see in the end of the 
present section (see Fig. [7] and its description) that this 
is indeed the case for our choice of initializations. Let us 
describe our initial state in detail. 

1. The initial space coordinates for the quarks, anti- 
quarks, and gluons are chosen at random within 
the finite box. 
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FIG. 4: (Color online) The reaction rates for elastic parton scattering (dashed green lines), gluon splitting (solid blue lines), and 
flavor-neutral qq fusion (short-dashed red lines) as functions of time for systems at different energy densities initially slightly 
out of equihbrium. (a) e = 1.1 GeV/fm^; (b) e = 4.72 GeV/fm^ 



2. The spectral properties (pole masses and the 
widths) of the quarks, antiquarks, and gluons are 
initially taken in the simple Lorentzian form (23) 
with two parameters for each parton type (M, T). 
Note that in the DQPM model one also assumes 
Lorentzian shapes for the parton spectral functions; 
however, we choose to start the system evolution 
not from the DQPM equilibrium spectral functions. 
For this purpose we deliberately employ an average 
value for the pole mass parameter in the spectral 
function of the strange quark at initialization (i.e., 
we choose Mg = Mu = Md). The other parameters 
[Mu, Md, Mg, Ti) are initially as in the DQPM. 
The spectral functions of the partons then evolve 
dynamically in time and in the final state may de- 
viate noticeably from the initial ones. We will see 
in the results of Section IIVI that indeed in the fi- 
nal thermalized state the dynamical gluon spectral 
functions deviate from the Lorentzian input and 
thus are not described by the DQPM ansatz. On 
the other hand, the pole mass of the strange quark 
dynamically reaches the correct value in equilib- 
rium. We stress here the importance of using off- 
shell transport for our studies. Only in case of the 
generalized transport propagation can we study the 
evolution of the spectral functions! 

3. We expect that in the chemically equilibrated state 
the ratio of strange quarks to the number of light 
{u or d) quarks is governed by the ratio of their 
masses (their flavor decomposition). We start our 
simulation from a flavor ratio, which is far from 
equilibrium; i.e, in the initial state the ratio of the 
number of s quarks to the number of u quarks and 
to the number of d quarks as 1:3:3 is taken as such 
that the strangeness is clearly undersaturated ini- 
tially. 



4. The initial momentum distributions and abun- 
dances of partons are given by the thermal distri- 
butions 



f{uj,p) = Cip'^ujpi{uj,p)np(^B)i^/Tinit) 



(58) 



where pi(a;,p) are the spectral functions (with i — 
q, q, g) and ngi^p^{uj /Tinit) are the Bose (Fermi) dis- 
tributions with a "temperature" parameter Tinu, 
which should not be misidentified with the final 
temperature T, which will be characteristic for the 
energy distributions of the particles after the ther- 
malization. The latter, "true" temperature T is 
well defined for the final, thermalized state, and in 
Sec. IIVI it will be extracted from the final particle 
spectra by fitting their high-energy tails. We will 
use this extracted final temperature T to study the 
equation of state of the partonic matter in PHSD 
in Sec. IIVI On the other hand, the value of the 
"temperature" parameter Tinu of the initial energy- 
momentum distributions and the numbers of par- 
tons (determined by the coefficients Ci) just define 
the total energy of the system (and in equilibrium 
the quark chemical potentials). 

5. The dynamical quarks, antiquarks, and gluons 
within the PHSD interact also via the mean fields. 
Note that the potential energy of this interaction is 
taken into account at initialization, so that it con- 
tributes to the total energy density. The strength 
of the quark and gluon potential energy in PHSD is 
given by the spacelike part of the 00 components of 
the energ y rn omentum tensor r'^° as in the DQPM 
(see Ref. '^). 

In the course of the subsequent transport evolution 
of the system by PHSD, the numbers of gluons, quarks. 
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FIG. 5: (Color online) Abundances of the u (solid red lines), d (short-dashed black lines), and s (dash-dotted blue lines) quarks 
-|- antiquarks and gluons (dashed green lines) as functions of time for systems at different energy densities, (a) £ — 1.1 GeV/fm''; 
(b) e = 4.72 GeV/fm^ 



and antiquarks change dynamically through inelastic and 
elastic collisions to equilibrium values. We observe in 
Fig. H that after about 20 fm/c (for e = 1.1 GcV/fm^) 
or 3 fm/c (for e = 4.72 GeV/fm^) the reactions rates 
are practically constant and obey a detailed balance for 
gluon splitting and qq fusion. In Fig. 2] the reaction rates 
for elastic parton scattering (dashed green lines), gluon 
splitting (solid blue lines), and flavor- neutral qq fusion 
(short-dashed red lines) are presented as functions of time 
at energy densities of 1.1 and 4.72 GeV/fm"^. We find 
that the rate of inelastic collisions relative to the elastic 
rate is larger at higher energy density; this is due to a 
larger gluon fraction with increasing energy density (or 
temperature) since gluons are more suppressed at low 
temperature due to their larger mass difference relative to 
the quarks. Nevertheless, it is worth mentioning that the 
elastic scattering between partons dominates in PHSD. 

A sign for chemical equilibration is the stabilization of 
the numbers of partons of the different species in time 
for i — >■ (X). In Fig. [S] we show the particle abundances of 
u, d, and s quarks-|-antiquarks (solid red, short-dashed 
black, and dash-dotted blue lines, respectively) and glu- 
ons (dashed green lines) for systems at energy densities 
of 1.1 and 4.72 GeV/fm^, which are above the critical en- 
ergy density (as in the previous figure) . We note in pass- 
ing that energy conservation within PHSD holds with an 
accuracy better than 10"'^ in these cases, which is a nec- 
essary requirement for our study. The slow increase of the 
total number of strange quarks and antiquarks during the 
time evolution reflects long equilibration times through 
inelastic processes involving strange partons. These time 
scales are significantly larger than typical reaction times 
of nucleus-nucleus collisions at SPS or RHIC energies. 
Note, however, that the rapidity and transverse momen- 
tum spectra of strange hadrons are well described by 
PHSD from lower SPS to top RHIC energies [Uli^. 
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FIG. 6: (Color online) The reactions rates for gluon splitting 
to pairs of strange quarks and antiquarks (solid blue line) and 
flavor-neutral ss fusion (short-dashed red line) as functions of 
time for a system at an energy density of 1.1 GeV/fm^. 



These findings appear to be in contradiction; however, 
the time scales from the box calculations cannot directly 
be applied to nucleus- nucleus collisions since the initial 
conditions are very different. The initial state in the box 
is chosen close to thermal parton equilibrium. This sup- 
presses the production of strange quark-antiquark pairs 
due to kinematics or available energy. The strangeness 
production in A -I- A collisions occurs mainly in the early 
stage oi A + A reactions where the system is rather far 
away from local thermal equilibrium and kinematical (en- 
ergy) constraints are subleading, i.e., particle collisions 
with large center-of-mass energies take place. These en- 
ergies are much larger than those in local thermal equi- 
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FIG. 7: (Color online) Abundances of the u, s quarks and 
gluons as functions of time for systems at e — 1.48 GeV/fm'^ 
with the different initial flavor ratios u: d: s — 3:3: 1 (solid 
lines), u:d:s = 2.5:2.5: 1 (dash-dotted lines), u:d:s = 2:2:1 
(dashed lines) and u:d:s = 1:1:0 (short-dashed lines). 



on the initial flavor ratios. This implies that our calcu- 
lations are stable with respect to the different initializa- 
tions, confirming that the system does reach equilibrium 
in our microscopic PHSD calculations. Since the equi- 
librium state is well defined by the PHSD calculations 
at each energy density (e.g., for times t > 120 fm/c), 
we may now proceed to study further properties of the 
system in dynamical equilibrium. 



IV. PHSD EQUILIBRIUM CALCULATIONS IN 
COMPARISON TO THE DQPM 

To compare the particle properties in the equilibrated 
dynamical model and in the DQPM, which has been de- 
veloped to describe QCD in equilibrium, we calculate dy- 
namically the different parton spectral functions. Let us 
consider the scalar parton density function ps defined (in 
equilibrium) by 



PsiT) ^ dg —- 



^ / 2tt J {2tt)^'^ 



,2vVPs(^,p)"s(c^/r)e(p2) 



librium, which makes the strangeness production more 
effective in A -I- A collisions and leads to lower strangeness 
equilibration times. Note that these arguments are sup- 
ported by the calculations also in HSD, where both col- 
liding and produced particles are hadrons [which happens 
at Alternating Gradient Synchrotron (AGS) and lowest 
SPS energies], as well as in PHSD, where the degrees of 
freedom are quarks and gluons [at top SPS, RHIC, and 
Large Hadron Collider (LHC) energies]. 

In Fig. [6] we present the time evolution of the reaction 
rates for gluon splitting to pairs of strange quarks and 
antiquarks (solid blue line) and ss fusion (short-dashed 
red line) for a system at an energy density of 1.1 GeV/fm"^ 
with the s and s quarks initially suppressed by a factor 
of 3 with respect to the light quarks. Accordingly, the 
initial rate for s-l-s — > g is suppressed by about a factor of 
9 and a large time for chemical equilibration is observed 
again. 

The results, shown in Figs. HHSl correspond to the ini- 
tial ratios between u, d, and s quarks (and antiquarks) 
taken as 

u:d: s — 3:3:1. 

We now vary the initial flavor decomposition and see if 
the system approaches the same final state (at constant 
energy density). In Fig. [7] we show the particle abun- 
dances of the M, and s quarks and gluons as functions of 
time for systems populated with the different initial flavor 
ratios: u:d:s — 3:3:1 (solid lines), u:d: s — 2.5 : 2.5 : 1 
(dash-dotted lines), u:d:s — 2:2:1 (dashed lines), and 
u\ d: s = 1:1:0 (short-dashed lines) while preserving 
the same energy density of the system e = 1.48 GeV/fm^ 
in all cases. One can see that the equilibrium values of 
the parton numbers for different flavors do not depend 



+ d 



'Qig) 



^y (^2v^p,(g-)(^,p) (59) 



X {npiiuj - Pg)/T] + UFiiLu + Ai,)/T]}e(p2), 

where ub and np denote the Bose and Fermi distribu- 
tions, respectively, while fj,q stands for the quark chemical 
potential. Here the scalar parton density is summed over 
gluons, quarks, and antiquarks. The number of gluonic 
degrees of freedom is dg = 16, while the fermion degrees 
of freedom amount to da 



iq = 2NcNf = 18 in case 



of three flavors {Nf = 3). The function 8(p ) (with 
p^ = w^ — p^) projects on timelikc four- momenta since 
only this fraction of the four-momentum distribution can 
be propagated within the light cone. In Eq. (|59p the par- 
ton spectral functions pj (with j = q, q, g) are no longer 5 
functions in invariant mass squared but taken as in (j30l) . 
Then the total number of timelike gluons g (quarks q or 
antiquarks q) in equilibrium (for p,q = 0) is given by the 
vector densities in thermodynamic equilibrium multiplied 
by the volume V: 
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(60) 



Note that for the scalar densities the integrand is the 
invariant mass divided by the energy to {y^/uj), while 
for the vector densities the integrand is simply 1 . For the 
energy spectrum we have 
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FIG. 8: (Color online) The spectra of the u and s quarks and gluons in equilibrium for different energy densities from the 
PHSD simulations (solid red lines) in comparison to the DQPM model (dashed blue lines). 



By choosing the momenta of the partons in the (narrow) 
interval |p| G [p_,p+], the energy spectrum is given by 
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^ ^Pg{q,q)i^^Prmd)nB(F)i'^/T), 



where pmid = ip+ — P-)/2 is the average momentum in 
the bin. 

In the transport approach we can construct the dis- 
tribution of partons with given energy and momentum 
as 
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which can be easily evaluated within the PHSD simu- 
lations in the box. Its counterpart within the DQPM 



model is 
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In Fig. m we show d^N/diudp for u and s quarks and 
gluons obtained by the PHSD simulations (red solid lines) 
of "infinite" partonic systems at energy densities of 1.1 
and 4.72 GeV/fm^. For comparison, we present on the 
same plots the DQPM assumptions (dashed blue lines) 
for the respective distributions. One can see that the 
DQPM distributions are in good agreement with the dy- 
namical calculations within PHSD for all quarks but de- 
viate from the simulations at high energy density for glu- 
ons. We will return to this apparent deviation below. 

Due to the off-shell dynamics in PHSD (cf. Sec. II.A) 
we have also access to the dynamical spectral functions in 
and out of equilibrium. Here we focus on the equilibrium 
state. Accordingly, we can compare the spectral func- 
tions of partons within the PHSD simulations in the box 
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FIG. 9: (Color online) The spectral functions of the u and s quarks and gluons in equilibrium from the PHSD simulations 
(solid red lines) for different energy densities in comparison to the DQPM model (dashed blue lines). 



and with the DQPM assumption for the spectral func- 
tions (j30p . Using the expression for the energy spectrum 
dMI, we get 
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In Fig. |9] we show the dynamical spectral functions 
p{u}) for u and s quarks and gluons as obtained by the 
PHSD simulations (red solid lines) for "infinite" partonic 
systems — at energy densities of 1.1 and 4.72 GeV/fm'^ — 
and the DQPM assumptions (dashed blue lines) for the 
spectral functions (|30l) at the corresponding energy den- 
sities of the system. 

We find that the dynamical spectral functions of 
quarks and gluons are generally fairly well described by 
the DQPM form ([50]) . However, there is a slight deviation 



visible at high energy density, especially for gluons. This 
deviation explains the difference between the dynamical 
results and the DQPM in Fig. \8\ The origin of the devi- 
ation can be traced back to the inelastic collisions of qq 
pairs forming gluons (P5|) in dense systems. The reactions 
favor the high-mass part of the gluon spectral function 
and predominantly populate dynamically the right-hand 
side from the gluon pole mass since the sum of the pole 
masses of quarks and antiquarks is larger than the pole 
mass of the gluon [cf. (|29l) ]. Indeed, let us recall that the 
inelastic collisions are more important at higher energy 
densities (cf. Fig.|4|). Moreover, from Fig.jTOlwe see that 
the elastic scattering rate of gluons is lower than that of 
quarks. Therefore, the inelastic interaction contributes 
considerably to the shape of the spectral function of glu- 
ons at high energy density, while it is not so important 
for the quarks at e = 1.1 GeV/fm^. In the DQPM it 
is assumed that the width in the spectral function is in- 
dependent of the mass, which indeed is found to be a 
good approximation if elastic scatterings dominate (as in 
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FIG. 10: (Color online) The total reaction rate for parton 
elastic scattering (short-dashed burgundy line) and separately 
the reaction rates for gluon-gluon (solid red line), gluon- 
quark (dashed blue line), quark-gluon (dotted green line), and 
quark-quark (dash-dotted dark yellow line) elastic scatterings 
as functions of time for the system at e = 1.1 GeV/fm'^. 
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FIG. 11: (Color online) The spectrum of u quarks in equilib- 
rium obtained by the PHSD simulations (red solid line) for 
systems at energy density 4.72 GeV/fm'^ in comparison to the 
thermal distributions with different temperatures: T — 243 
MeV (dashed blue line), T = 223 MeV (dash-dotted green 
line), and T = 263 MeV (short-dashed burgundy line). 



case of the quarks and antiquarks). However, the inelas- 
tic interaction of partons in PHSD is dominated by the 
resonant gluon formation, which dynamically generates 
a mass-dependent width for the gluon spectral function. 
This dynamical effect in the gluon width is not incorpo- 
rated in the DQPM assumption pop . Accordingly, the 
PHSD simulations for systems in equilibrium supersede 
the DQPM assumptions but well reproduce the DQPM 
assumptions in the fermionic sector. 

Note that the calculations oicPN/dtudp in PHSD in the 
box in the final, equilibrated state allows us to extract 
the temperature of the "infinite" parton matter. We ob- 
tain the final temperature by fitting the parton spectrum 
obtained by the PHSD simulations with the product of 
the Bose (Fermi) distribution and a (Lorentzian) spec- 
tral function [cf. formula (IM|) ] . In Fig. [Til we show the 
spectrum of u quarks from the PHSD simulations (solid 
red line) for a system at energy density 4.72 GeV/fm'^ in 
comparison to the fit with different temperatures: T = 
243 MeV (dashed blue line), T = 223 MeV (dash-dotted 
green line), and T — 263 MeV (short-dashed burgundy 
fine) . All three curves were normalized to coincide at the 
peak of the spectral function. One can see that the high- 
momentum behavior of the distribution is governed by 
the temperature and that the temperature T = 243 MeV 
gives the best fit at the energy density of 4.72 GeV/fm'^. 
We note that the same procedure is repeated for each 
particle species and each value of the energy density. 

The question of whether the equation of state from 
the PHSD in equilibrium compares reasonably with the 
lattice data from Ref. [33| can now be aswered. To this 
end we present in Fig. [T^] the equation of state extracted 
from the PHSD calculations in the box (red stars) in 
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FIG. 12: (Color online) The scaled energy density e/T* (red 
stars) and the potential energy fraction of the scaled energy 
density (green open squares) extracted from the PHSD calcu- 
lations in the box in comparison to the IQCD data from Ref. 
[331 (blue open circles). 



comparison to the respective results from the Wuppertal- 
Budapest group [33| (blue open circles) as functions of 
the temperature T. We also show the potential energy 
contribution to the equation of state extracted from the 
PHSD calculations in the box (green open squares) that 
is equivalent to the DQPM potential energy density. We 
find that the equation of state implemented in PHSD is 
well in agreement with the DQPM and the IQCD results. 
This finding implies that PHSD dynamically describes 
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FIG. 13: (Color online) Abundances of u, d, and s quarks (solid red lines, short-dashed black lines, and dash-dotted blue 
lines, respectively) and u, d, and s antiquarks (dotted red lines, short-dotted black lines, and short-dash-dotted blue lines, 
respectively) and gluons (dashed green lines) as functions of time for systems at e = 1.48 GeV/fm^ and at different quark 
chemical potentials: (a) fiq = 50 MeV; (b) Hq = 150 MeV. 
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FIG. 14: (Color online) The spectrum of u quarks in equilibrium as obtained by the PHSD simulations in the box (red solid 
lines) in comparison to the DQPM model (dashed blue lines) for systems at an energy density of 1.48 GeV/fm'^ at different 
quark chemical potentials: (a) fiq — 50 MeV; (b) ^q = 150 MeV. 



systems of quarks, antiquarks, and gluons in equilibrium 
that have the same properties as explicit QCD calcula- 
tions on the lattice. 



V. FINITE QUARK CHEMICAL POTENTIALS 

We have seen in the previous section that the dy- 
namical calculations within PHSD reproduce equilibrium 
properties of QCD matter as seen in lattice QCD cal- 
culations at fj,q = 0. Let us now proceed further and 
study within the dynamical approach the quark and 
gluon properties at finite quark chemical potential fiq, 
which are currently not yet well established in lattice 



QCD calculations. 

In Fig. [13] we present the particle abundances of u, d, 
and s quarks (solid red, short-dashed black, and dash- 
dotted blue lines, respectively), u, d, and s antiquarks 
(dotted red, short-dotted black, and short-dash-dotted 
blue lines, respectively) and gluons (dashed green lines) 
as functions of time for systems at e = 1.48 GeV/fm'^ and 
at quark chemical potentials fXq of 50 and 150 MeV. Again 
chemical equilibrium is achieved for large times but now 
the abundances of quarks and antiquarks differ consid- 
erable (especially for /Xg = 150 MeV). A closer inspec- 
tion also shows that the strangeness equilibration pro- 
ceeds slower since the amount of flavor-neutral qq pairs 
decreases with increasing fig. Note that the gluon abun- 
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dance in the equilibrium stage does not depend on the 
initiahzation. 

The phase boundary Tdfiq) in the DQPM (and PHSD) 
is defined by demanding that the phase transition hap- 
pens at the same critical energy density Sc for all fiq. 
The prediction of Fig. [13] might in future be compared to 
IQCD calculations at finite fiq once physical quark masses 
can be incorporated into the IQCD calculations. 

In Fig. [14] we show d'^N/dcodp for u quarks obtained 
by the PHSD simulations (red solid lines) of "infinite" 
partonic systems at e = 1.48 GeV/fm^ and at quark 
chemical potentials of 50 and 150 MeV. For comparison, 
we present on the same plots the DQPM assumptions 
(dashed blue lines) for the respective distributions. The 
agreement is fairly good since the inelastic channels are 
further suppressed with increasing fig. Note that in the 
present version the DQPM and PHSD treat the quark- 
hadron transition as a smooth crossover at all fiq. There 
are, however, some physical arguments in favor of a first- 
order phase transition at large /ig and for the existence 
of a critical endpoint for the first-order transition line in 
the T-^q plane. Presently, we are not able to calculate 
the properties of a quark-gluon system close to a criti- 
cal endpoint. It is also not yet clear whether such an 
endpoint exists. 



VI. SCALED VARIANCE, SKEWNESS AND 
KURTOSIS 

In this section we address higher moments of parton 
distributions in the sQGP within the PHSD approach 
and study the equilibration of fluctuation observables as 
well as the size of fluctuations in equilibrium. We recall 
that various fluctuation observables have been addressed 
theoretically within IQCD |42l - [45| as well as within effec- 
tive models |46l - l5l| . Furthermore, some of these observ- 
ables have been studied experimentally by the various 
collaborations at the SPS and at RHIC. Most of these 
have been evaluated in the HSD approach including the 
individual detector acceptance and experimental biases. 
For a recent review we refer the reader to Ref. [S^l (and 
references cited therein). The evaluation of the various 
fluctuations in PHSD is straightforward and in this sec- 
tion performed for ^q = 0. 



A. Scaled variances 

We start with scaled variance 

^2 



and cr^ is the sample variance given by 
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where ^ is the mean value of the observable x averaged 
over N events, 



1 ^ 
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(68) 



The standard error of the scaled variance uj is given by 
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and TO4 is the fourth central moment. 
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In Fig. [15] we show the scaled variances lo for par- 
ticle number fluctuations as functions of time for the 
quarks-|-antiquarks of all flavors (green open circles) , for 
separate quark flavors — u (red squares), d (blue down- 
triangles), and s (green up-triangles) — and for gluons 
(red opened squares), for a system at an energy den- 
sity of 1.48 GeV/fm'^. The same results are presented 
in Fig. [in]for a system at e = 4.72 GeV/fm^. Note 
that in the grand canonical ensemble, i.e., for an equi- 
librium system with constant temperature (due to the 
presence of a thermostat) and with thermal fluctuations 
of the total system energy, one would expect w « 1 for all 
particle number fluctuations. On the other hand, for an 
isolated statistical system the global energy conservation 
for the microcanonical ensemble leads to a suppression of 
the particle number fluctuations and thus to uj < 1 (see 
Ref. ^3] for more details). As seen from Figs. [T5l and [TBI 
the equilibrium values of to are smaller than 1 . This can 
be interpreted as a consequence of the total energy con- 
servation, which is flxed rather more strictly (but still nu- 
merically not exactly) in our PHSD box calculations than 
that in the grand canonical ensemble. In a mixture of dif- 
ferent particle species the influence of the global energy 
conservation on particle number fluctuations is different 
for different species in the mixture. The suppression ef- 
fects are stronger for those species that contain larger 
fractions of the total system energy. The scaled variance 
of all charged particles (i.e., quarks plus antiquarks of all 
flavors) is lower than that of gluons or of a single quark 
flavor. This reflects the larger energy fraction stored in 
all quarks. For illustration, we show in Fig. [T7] the total 
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FIG. 15: (Color online) The scaled variance as functions of time for a system at e = 1.48 GeV/fm for (a) all charged particles 
(green open circles) and giuons (red open squares) and (b) different quark flavors: u (red squares), d (blue down triangles), 
and s (green up triangles) quarks + antiquarks. 
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FIG. 16: (Color online) The scaled variance as functions of time for a system at e = 4.72 GeV/fm^ for (a) all charged particles 
(green open circles) and giuons (red open squares) and (b) different quark flavors: u (red squares), d (blue down triangles), 
and s (green up triangles) quarks + antiquarks. 



energy of partons (dash-dotted green lines), the energy 
of all charged partons (solid red lines) , and the energy of 
giuons (dashed blue lines) as functions of time for sys- 
tems at energy densities of 1.48 and 4.72 GeV/fm-^. We 
observe that the scaled variances reach a plateau in time 
for all observables and energy densities. The scaled fluc- 
tuations in the gluon number are more pronounced at 
e = 1.48 GeV/fm^ since the fraction of the gluon energy 
is quite small at this energy density. The difference with 
respect to the scaled variance of all charged partons de- 
creases with energy due to the higher relative fraction of 
the gluon energy, as discussed before. Due to the ini- 
tially lower abundance (thus smaller energy fraction) of 
strange quarks the respective scaled variance is initially 



larger but reaches the same asymptotic value as the light 
quarks in the course of the time evolution. Accordingly, 
the fluctuations in the fermion number are flavor blind 
in equilibrium. 

It is interesting to study the scaled variance for a cell 
inside the box as a function of the cell volume. This 
can easily be achieved by subdividing the total volume 
V — 9'^ fm'^ in different subvolumes Vn of equal size and 
evaluating the scaled variance in each subvolume. Fi- 
nally an average over the n subvolumes Vn is performed. 
In Fig. [18] we present the scaled variance as functions of 
n — V/Vn in the box for all charged particles (red open 
circles) and giuons (blue open squares) and for different 
quarks flavors [u (red open squares), d (blue open cir- 
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FIG. 17: (Color online) The total kinetic energy of partons (dash-dotted green lines), the kinetic energy of all charged partons 
(solid red lines), and the kinetic energy of gluons (dashed blue lines) as functions of time for systems at different energy 
densities: (a) e = 1.48 GeV/fm^'; (b) e = 4.72 GeV/fm^ 
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FIG. 18: (Color online) The scaled variances in equilibrium (at t > 120 fm/c) as functions of relative system size n — V/Vn 
at energy densities of 1.48 and 4.72 GeV/fm'^, where V is the default box volume and Vn is the subsystem volume for (a) all 
charged particles (red open circles) and gluons (blue open squares) and (b) different quark flavors: u (red open squares), d 
(blue open circles), and s (green open triangles) quarks + antiquarks. Note that n = 1 corresponds to a subsystem volume 
Vi = 1/ = 9^ = 729 fm^; n = 10 stands for Vio = 72.9 fm^; while n = 0.2 means a system of volume Vb.2 = 5 x 729 = 3645 fm^, 
which is larger than our default box size V. 



cles), and s (green open triangles) quarks + antiquarks] 
for systems at energy densities of 1.48 and 4.72 GeV/fm'^, 
respectively. The inserts show the observablcs for larger 
box sizes by up to about a factor of 8 (n « 0.15) in order 
to explore the thermodynamic limit. Indeed, our calcu- 
lations demonstrate that the scaled variances no longer 
change (within statistics) when increasing the volume of 
the box by up to about an order of magnitude, thus 



approaching the thermodynamic limit. We recall that 
cj = 1 for a Poissonian distribution (dash-dotted black 
lines). The impact of total energy conservation in the 
box volume V is relaxed in the subvolume Vn- This in- 
fluence becomes weaker for n ^ 1, i.e., Vn ^ V. There- 
fore, in the subvolume Vn the energy fluctuates and these 
fluctuations behave as in the grand canonical ensemble 
for n ^ 1. The remaining part of the box plays — in 
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FIG. 19: (Color online) The skewness as a function of time for a system at e = 1.48 GeV/fm^ for (a) all charged particles 
(green open circles) and giuons (red open squares) and (b) different quarks flavors: u (red squares), d (blue down triangles), 
and s (green up triangles) quarks. 



this limit — the role of a thermostat for the cell Vn ■ This 
explains the behavior a; = 1 for all scaled variances at 
large n as seen in Fig. [181 Such a behavior can be also 
expected from the "law of rare events" : the scaled vari- 
ances for all observables approach the Poissonian limit 
when one considers only a tiny fraction of all particles in 
the system. 

This observation raises a new question concerning the 
event-by-event fluctuations in nucleus-nucleus collisions 
within a viscous hydrodynamical approach. The basic 
requirement of this approach is that the local cell size — 
in which a possibly chemical and kinetic equilibrium is 
achieved — is small compared to the macroscopic dimen- 
sion of the system; in particular, the gradients in the 
energy density should be small. In each cell then equilib- 
rium values for averages as well as fluctuations of observ- 
ables should be considered within the grand canonical 
treatment. Thus, the influence of the conservation laws 
(both energy-momentum and charge conservation) gets 
lost. However, the influence of the global conservation 
laws on fluctuation observables is by no means negligible 
even in the thermodynamical limit, if the detector would 
accept an essential fraction of all particles. 



B. Skewness 

The skewness [SJ] characterizes the asymmetry of the 
distribution function with respect to its average value. 
If the bulk of the data are at the left and the right tail 
is stretched out, then the distribution is skewed right or 
positively skewed; if the peak is toward the right and 
the left tail is more pronounced, then the distribution 
is skewed left or negatively skewed. The definition of 



skewness is as follows: 
51 = 
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where TO2 and rrij, are the second (variance) and third 
central moments, respectively, with 
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The skewness of a sample is given by 
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In Fig. [T9| we show the skewness as functions of time 
for all charged particles (green open circles) and giuons 
(red open squares) and for different quarks flavors \u (red 
squares), d (blue down triangles), and s (green up trian- 
gles) quarks] for a system at e = 1.48 GeV/fm'^. Note 
that the skewness is equal to zero for symmetric distri- 
butions (dash-dotted black lines). We find that in our 
case the skewness of the number of all charged particles 
tends to be antisymmetric to the skewness of the num- 
ber of giuons, but both are compatible with zero for the 
present accuracy of the calculations. We only show the 
results for a single energy density since our findings are 
independent of the energy density. 



C. Kurtosis 

The height and sharpness of the distribution peak rel- 
ative to a number is called kurtosis ^54)] . Higher values of 
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FIG. 20: (Color online) The excess kurtosis as a function of time for a system at an energy density of 1.48 GeV/fm^ and the 
corresponding IQCD results (dashed orange lines) from Ref. [4J] for (a) all charged particles (green open circles) and giuons 
(red open squares) and (b) different quarks flavors: u (red squares), d (blue down triangles), and s (green up triangles) quarks. 



kurtosis indicate a higher, sharper peak; lower values of 
kurtosis indicate a lower, less distinct peak. The kurtosis 
is defined as 



_ m4 _ m4 
P2 — ^ - ^ 
1712 cr 



(77) 



where 1714 is determined by (I72p . It is equal to 3 for a 
normal distribution, so often the excess kurtosis is pre- 
sented which characterizes the deviation from a normal 
distribution, 



The sample excess kurtosis then is defined by 
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The standard error of the kurtosis is given by 
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where AGi is determined by (1761) . 

In Fig. [20] we present the excess kurtosis as functions 
of time for all charged particles (green open circles) and 
giuons (red open squares) and for different quarks flavors 
[u (red squares), d (blue down triangles), and s (green 
up triangles) quarks] for a system at an energy density of 
1.48 GeV/fm'^. Note that the excess kurtosis is equal to 
zero for normal distributions (dash-dotted black lines). 
The IQCD results from Ref. [4J| are nonzero and shown 
by the dashed orange lines. We find that in our case 
the excess kurtosis of the number of all charged particles 
is equal to the excess kurtosis of the number of giuons. 
However, within statistical errors, the excess kurtosis is 
compatible with zero as well as with the IQCD results 
from Ref. [4J| for giuons and charged particles. This 



VII. EQUILIBRATION TIMES 

An inspection of the time evolution of the scaled vari- 
ances in Figs, [in] and [THJ shows that the equilibration of 
the various scaled variances occurs on time scales that 
are shorter than the time scales for the equilibration of 
the average values of the observables. In order to quan- 
tify this observation we fit the explicit time dependence 
of the abundances and scaled variances by the function 



Ot = Ot=o + [Ot^oo ~ Ot^o) (1 - e-*/"-) 



(81) 



finding holds for all energy densities considered. 



which defines a characteristic equilibration time T^q . The 
results of our fits for different observables and energy den- 
sities are summarized in Table |I] For all particle species 
and energy densities, the equilibration time r^q is found 
to be shorter for the scaled variances than for the aver- 
age values. This is most pronounced when considering all 
charged partons but less distinct for strange quarks. Ac- 
cordingly, scaled variances may achieve an equilibrium, 
even if the average values of an observable are still out 
of equilibrium. This finding is reminiscent of strongly 
interacting quantum systems evaluated on the basis of 
Kadanoff-Baym equations in Ref. [24,], where quantum 
fluctuations stabilize early in time, i.e., long before a ki- 
netic or chemical equilibrium is achieved. 



VIII. SUMMARY AND CONCLUSIONS 

In this work we have employed the PHSD off-shell 
transport approach to study partonic systems slightly 
out of equilibrium as well as in equilibrium in a finite 
box with periodic boundary conditions, thus simulating 
"infinite" partonic matter. After a brief recapitulation 
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Particle type 


Equilibration times Teq fm/c 


£ = 1.48 GeV/fm=* 


e = 4.72 GeV/fm» 


abundance 


UJ 


abundance 


CJ 


u + u 


43 


16 


21 


6 


d + d 


45 


14 


21 


5 


s + s 


35 


28 


19 


17 


gluons 


18 


5 


18 


4 


all charged 


18 


3 


18 


2 



TABLE I: Equilibration times for the abundances and the 
scaled variances for the different particle species and two val- 
ues of the energy density. 



of off-shell dynamics in phase space we have specified in 
more detail the ingredients of PHSD, i.e., the retarded 
self-energies of the partons and the elastic and inelastic 
cross sections for partons. Furthermore, we have recapit- 
ulated the basic equation for the transition from partons 
to hadrons, i.e., the dynamical hadronization that incor- 
porates all conservation laws as well as an increase in 
total entropy for a rapidly expanding systems. 

We have demonstrated explicitly that partonic systems 
at energy densities e, above the critical energy density 
£c ~ 0.5 Gev/fm'^, achieve kinetic and chemical equi- 
librium in time. Furthermore, the energy density of the 
partonic system at fixed temperature and quark chemical 
potential for u^, = is well in line with the lattice QCD 
calculations |39| in equilibrium. This allows us to study 
explicit equilibration times for different observables when 
initializing the partonic system slightly out of equilibrium 
and also at finite but moderate quark chemical potential 
/ig ^ 0. Most strikingly, we find that the strangeness de- 
gree of freedom equilibrates on time scales that are large 
compared to the reaction times in relativistic nucleus- 
nucleus collisions. Nevertheless, the application of PHSD 
to these reactions from low SPS to top RHIC energies 
provides a good description of strangeness observables in 
rapidity and transverse momentum [2l|,|40l. At first sight 
this might look like a contradiction; however, the initial 
nucleon-nucleon collisions in relativistic A + A reactions 
occur at much larger invariant energies then those in local 
thermal equilibrium. This makes the strangeness produc- 
tion more effective in A+A collisions. Moreover, the kaon 
to pion ratio is enhanced at midrapidity — where statisti- 
cal model fits are performed — due to a narrower rapidity 
distribution of kaons relative to pions. This effect is seen 
experimentally and is also present in the nonequilibrium 
transport (HSD and PHSD) calculations. 

In addition to equilibration phenomena of average val- 
ues for observables such as particle number or charged 
particle number we have studied the dynamics of fluctu- 
ation observables in and out of equilibrium. For all ob- 
servables the equilibration time Teq is found to be shorter 
for the scaled variances than for the average values. This 



is most pronounced when considering all charged partons 
but less distinct for strange quarks. Accordingly, scaled 
variances may achieve an equilibrium even if the average 
values of an observable are still out of equilibrium. This 
finding is reminiscent of strongly interacting quantum 
systems evaluated on the basis of Kadanoff-Baym equa- 
tions 2J], where quantum fiuctuations stabilize early in 
time, i.e., long before a kinetic or chemical equilibrium is 
achieved. 

The scaled variances for the fluctuations in the num- 
bers of different partons in the box show an influence 
of total energy conservation. We observe a suppression 
of the parton number fluctuations in comparison to the 
fluctuations in the grand canonical ensemble. Further- 
more, by dividing the box into several cells we have cal- 
culated the scaled variances of different observables in 
the cell as functions of the cell size. The scaled vari- 
ances for all observables approach the Poissonian limit 
with a; = 1 when the cell volume is much smaller than 
that of the box. This observation indicates that global 
conservation laws (for energy-momentum and charges) 
are not important when one detects only a small frac- 
tion from all particles in the system. However, if the 
fraction of the accepted particles is comparable to that 
in the whole system, the influence of global conservation 
laws on fluctuation observables is not negligible even in 
the thermodynamic limit. We have shown, furthermore, 
that the scaled variances no longer depend on the size of 
the box when increasing it up to about an order of magni- 
tude up to ^5000 fm"^. Accordingly, the continuum limit 
has approximately been reached in the calculations. 

Our analysis of the skewness and kurtosis gives prac- 
tically vanishing values for these observables in time 
and especially in equilibrium within the limited statistics 
achieved. We mention that our results within statistics 
are also compatible with the IQCD results from Ref. [44| . 
This issue will have to be readdressed in future along with 
an evaluation of transport coefficients such as the shear 
viscosity and bulk viscosity as functions of temperature 
and quark chemical potential 55']. 
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